NASA TM X- 52761 



NASA TECHNICAL 

MEMORANDUM 


l 


«* 


£ 

C\J 

ir\ 

i 

X 



< 

CO 

< 



EXPERIMENTAL VERIFICATION OF EFFECTS OF TURBULENT 
BOUNDARY LAYERS ON CHEMICAL-KINETIC 
MEASUREMENTS IN A SHOCK TUBE 

by F. E. Belles and T- A. Brabbs 
Lewis Research Center 
Cleveland, Ohio 

1 

TECHNICAL PAPER proposed for presentation at 

J# Symposium on Combustion sponsored 

by the Combustion Institute 

Salt Lake City, Utah, August 23-29, 1970 



EXPERIMENTAL VERIFICATION OF EFFECTS OF TURBULENT 

BOUNDARY LAYERS ON CHEMICAL -KINETIC 

MEASUREMENTS IN A SHOCK TUBE 

by F. E. Belles and T. A. Brabbs 

Lewis Research Center 
Cleveland, Ohio 


TECHNICAL PAPER proposed for presentation at 

13 m Symposium on Combustion sponsored 
by the Combustion Institute 
Salt Lake City, Utah, August 23-29, 1970 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 



E-4938 


EXPERIMENTAL VERIFICATION OF EFFECTS OF TURBULENT BOUNDARY LAYERS 


ON CHEMICAL-KINETIC MEASUREMENTS IN A SHOCK TUBE 
by F. E. Belles and T. A. Brabbs 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

SUMMARY 

Experiments were conducted to demonstrate the existence and magnitude 
of the effects of nonuniform shock- tube flow on kinetic results. Atomic- 
oxygen concentration was observed in a rich H 2 / 02 /C 0 /Ar mixture behind 
incident shocks that produced turbulent boundary- layers and essentially con- 
stant shock-front temperatures, at eight initial pressures from 50 to 120 
torr. The exponential time-constant for growth of [o3 and the time at which 
peak concentration occurs are governed by the rate constant for 
H + O 2 — * OH + 0. Values of the rate constant were fit to the 1 20- torr 
data using two analyses: one incorporating the boundary- layer effects on 
flow properties and residence time as formulated by Mirels; and the other 
making the conventional assumption of constant post-shock conditions. Two 
sets of predictions were then made for the lower-pressure runs and compared 
with experiment. The comparison showed that the conventional analysis is 
incorrect; that there is a large error in the rate constant derived from it; 
and that Mirels' formulas satisfactorily account for the observed results. 

INTRODUCTION 

Mirels has suggested that flow non-uniformities induced by the boundary 
layer behind an incident shock may significantly affect shock-tube studies 
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of chemical kinetics. His discussion* of these effects was for the quasi- 
steady case in which the shock and the contact surface have reached their 
maximum separations^. Under these circumstances the temperature, density, 
and gas residence time are all predicted to increase with distance behind 
the shock front, even though attenuation of the shock is no longer considered. 

Mirels' suggestion has apparently not had a great impact on chemists 
using shock tubes. There have been only a few published instances in which 
corrections for changing conditions were mentioned and still fewer in which 
corrections were applied. And even in the most thorough analysis** that has 
been performed, there were no concurrent experiments to show that the kinetic 
data actually required correction, nor that Mirels' formulation* of boundary- 
layer effects was a proper basis for making corrections. Thus, it is not 
surprising that others should hesitate to revise the usual practice of inter- 
preting shock-tube data as if the flow remained uniform. 

The present work was done to show that the changing conditions behind 
shock waves do indeed have a significant effect on kinetic data, and to 
investigate how well Mirels' description of the flow applies to experiments 
done in the presence of turbulent boundary layers. The approach was to do 
a kinetic experiment in reverse: Rather than measure rate data, correct 
them for boundary^- layer effects, and then compare the results with other 
data which are themselves of uncertain accuracy, we exploited the sensitivity 
to changing flow conditions of a simple chemical system that is already well 
understood. A very rich H2/O2/CO mixture, when shocked to a suitable temper- 
ature, emits a brief pulse of blue light. The shape and time of occurrence 
of this pulse are governed almost entirely by a single reaction, the rate 
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and progress of which are sensitive to changes in temperature, density, and 
residence time. Therefore, the light pulse serves as a chemical transducer 
with which to probe the flow properties. 

A series of observed light pulses was analyzed with and without con- 
sideration of boundary- layer effects. The simplicity of the chemistry made 
it possible to carry out the two analyses by using in each a different value 
of the crucial parameter: the rate constant of the dominant reaction. The 
validity of each method of analysis could then be judged by its success in 
accounting for the experimental results. 

BACKGROUND 

Boundary-Layer Effects 

After the shock and contact surface have reached their maximum sepa- 
ration, the temperature, density, and residence time of the gas all increase 
behind the shock as they would in a subsonic flow if the cross-sectional 
area of the passage were increasing. For the case of a turbulent boundary 
layer, the effective area is given by the following equation' 

(Ax/A) = [l - (x/l m )0-8] (1) 

where 

A v effective cross-sectional area of tube 

A 

A actual cross-sectional area of tube 

x distance behind shock, equal to (laboratory time) x (shock velocity) 
l m maximum separation of shock apd contact surface. 

The value of l m depends on the gas mixture and its initial temperature and 
pressure, the hydraulic diameter of the tube, and the Mach number of the 
shock, and can be calculated from formulas given in reference 3* 



-4- 

Chemistry 

The Intensity of the blue light (CO flame-band radiation) that is 
emitted when a H 2 /O 2 /CO mixture is shock-heated to sufficiently high 
temperature is directly proportional to the product of CO and 0 concen- 
trationsS. If the mixture contains relatively large amounts of H 2 and CO 
compared to 02> very little of the CO is consumed; thus, the light intensity 
essentially varies with the 0 concentration, [oj . 

Wei 1-conf i rmed theory^"® shows that [oj is governed by the following 
sequence of reactions. Immediately after the gas has been shock-heated, 
one or both of the following initiation reactions produce a small concen- 
tration of chain-carriers (H, 0, or OH): 


H 2 + 02 

— » 20H 

(i) 

CO + O 2 

— > co 2 + 0 

(ii) 


The brief initiation period is followed by an induction period during which 
the concentrations of all three chain carriers grow exponentially with time 
and all with the same exponential time constant, of , as a result of the 
chain-branching reactions of hydrogen combustion^. Various of these reactions 
assume different degrees of importance in different mixtures but in very rich 
ones such as that used in the present work (H2/02= 1 0/ l ) , ©f is governed 
almost entirely by reaction (ll) of the H 2 - 0£ combustion mechanism: 

H + O 2 — * OH + 0 (II) 

In the limit of extreme H 2 /O 2 ratios and in the absence of boundary- layer 
effects, ©£ is given by the following equation^: 

e<r= 2f>2] k 2 [°2] 


( 2 ) 
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where 

cC exponential time constant for growth of 0 concentration in laboratory 
time-coordinate system 
p 2 j density ratio across the shock 
k 2 rate constant for reaction (l l) 

[O 2 I oxygen concentration 

The short supply of oxygen in very rich mixtures eventually causes the 
exponential buildup of ro] to be offset by rapid depletion®. As a result, 
flame-band radiation is observed to peak and then to drop rapidly. The 
time at which the peak occurs, 'V , is also governed mainly by the rate of 
reection (II) in the limit of extreme H 2 /O 2 ratios, although there is a minor 
effect of the rate of initiation by reaction (i) as well. This can be seen 
from the following equation, which is derivable from relations given in 
references 6 and 8: 

t 

where 

'C time at which peak 0 concentration occurs in laboratory time - coordinate 
system 

kj rate constant for reaction (i) 

C expression containing factors that remain constant in this experiment. 
Despite the small proportion of O 2 present in the extremely rich mixtures 
under consideration, both the exponential rise and the peak in [0] occur 
before any appreciable amount of it has been consumed. 
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APPROACH 

Experimental Conditions 

The experimental objective was to obtain a set of measurements of«C and 
X that would, when analyzed with and without consideration of effects 
induced by a turbulent boundary layer, reveal the presence of such effects. 
In order to accomplish this, several requirements had to be considered in 
planning the experiments. These requirements are listed below, together 
with the measures taken to fulfill them. 

Simple chemistry : To avoid the ambiguities that arise when experiments 

must be analyzed by permuting the rate constants of many reactions. 

Although Eqs. (2) and (3) are not quantitatively correct for a H2/O2 = 10/1 
mixture, they do correctly show that k2 is the only adjustable parameter of 
any consequence. 

Low oxygen concentration ; To prevent excessive heat release that might 
cause irregular waves^. The mixture chosen was H2/02/C0/Ar = 5/0.5/6/88.5. 
Not the slightest evidence of irregular waves was found, either in light- 
emission or pressure records. 

Low temperature : To keep k2 low and thereby delay the light pulse so 

that the boundary layer could exert significant effects (see Eq. (1)). 
However, the temperature must not be too low lest chain-breaking by 
H + O2 + M — * HO2 + M assume undue importance. 

Constant temperature : To make pressure the only experimental variable 

and to eliminate arguments about the influence of the activation energy 
chosen for (II). By tr ia 1 -and-error changes in driver pressure, shock-front 
temperatures were held to the narrow and suitably low range of 1 140° + 12°K. 
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High pressure : To assure turbulent boundary layers^. The lowest initial 

pressure used was £0 torr. Runs were made at 10-torr intervals up to 120 torr. 

Analysis 

After the oscilloscope records of light intensity had been plotted semi- 
logarithmically against laboratory time to yield values ofoC and "t , the 
data were analyzed by means of a computer program. This program integrated 
the equations of chemical change for reactions occurring behind a shock 
wave and was equipped to handle two cases: (1) the case in which the apparent 
area of the tube varies under the influence of the boundary layer in accor- 
dance with Eq. (1), with resulting changes in temperature, density, and 
residence time of the gas; and (2) the conventionally-assumed case in which 
the area remains constant and the only changes in properties are those due 
to chemical reaction. The two approaches will be referred to subsequently 
as the varying-area and the constant-area methods of analysis, respectively. 

The procedure was as follows. First, the run at the highest initial 
pressure (120 torr) was analyzed by both methods. By trial-and-error cal- 
culations, two different values of were found that would closely reproduce 
the observed «C . Although the chemistry is dominated by reaction II, 12 
other initiation ((i) and (ii)), chain-branching, and recombination reactions 
pertinent to the H 2 -O 2 -CO system were included in the analysis. Values for 
their rate constants were taken from recent literature and were identical 
for both methods of calculation. 

Next, the calculated peak times were matched to theT observed in the 
120-torr run by modifying kj (see Eq. (3)). This resulted in two different 
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values of kj corresponding to the two methods of analysis. 

The rate constants obtained in this way were as follows: 


k£ (constant area) 

= 2.10xl0 1/f exp(-l6600/RT) 

(4a) 

k 2 (varying area) 

= 1.44xl0 ,if exp(-l6600/RT) 

(4b) 

k| (constant area) 

= 2. 10xl0 ,2 exp(-39000/RT) 

(5a) 

K- (varying area) 

= 1.20xl0 12 exp(-39000/RT) 

(5b) 


Only the pre-exponential parts of Eqs. (4) and (5) resulted from the trial- 
and-error fits; the activation energies were assigned. That for reaction 
(ll) was derived from an Arrhenius plot, covering a wide temperature range, 
of rate constants from many literature sources that were considered reliable. 
The activation energy assigned to reaction (i) is the value reported in 
reference 10. 

The final step in the analysis was to use the two sets of rate constants 
to predict constant- and varying-area values of oC and 'C for the runs made at 
the pressures below 120 torr. Comparison of the predictions with the experi- 
mental data could then be made. 

APPARATUS AND PROCEDURE 
Shock Tube 

The tube was a single piece of stainless steel, 5.7 meters long. The 
internal dimensions were 6.4 x 6.4 cm. with corners rounded to a radius of 
1.3 cm. The entire length of the tube was ground to constant inside dimen- 
sions and then honed to a highly-polished finish. 

Stations for shock-wave detectors were located at 15-cm, intervals in 
the downstream portion of the tube. A piezoelectric pickup which triggered 
a raster oscilloscope was followed by four matched pressure transducers 
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for velocity measurements. These transducers were of the quartz piezoelectric 
type and had short rise times. 

Midway between the last two stations was a pair of 2.5-cm. diameter 
windows made of calcium fluoride and located opposite one another. A thin- 
film gauge was located to provide an accurate indication of the time at 
which a shock wave arrived at the center of the windows and to ascertain that 
the boundary layer was turbulent. All pickups and windows were carefully 
installed with their surfaces flush with the inner wails of the tube. 

The assembled tube could be evacuated to a pressure of about 1 micron 
and had a leak rate less than 0.2 micron/minute. A liquid nitrogen cold 
trap in the vacuum line guarded against the possible back-migration of pump 
oil. 

Velocity Measurement 

The overriding experimental requirement was the precise measurement of 
shock velocity. Uncertainties as large as 1 percent, frequently tolerated 
in shock-tube work, would have completely vitiated the results of this 
particular experiment. This can be seen by examining Fig. 1. These com- 
puted profiles of the product of CO and 0 concentrations are equivalent to 
light-emission histories. The differences in shape and peak time are 
readi ly apparent. 

The required precision was obtained by electronically processing the 
signals from the four matched transducers so as to produce pulses. These 
were displayed on a raster oscilloscope along with 1-microsecond timing 
marks. In this way, velocities were measured with an uncertainty of 0.2 


percent. 
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Light-Detection System 

Flame-band emission was observed through one of the windows by means 
of a 1:1 optical transfer system. A lens imaged the center of the shock 
tube onto a slit, 0.25 mm wide, which acted as a field stop. The illuminated 
slit became the light source for a spherical mirror, which transferred the 

o 

luminous image at unit magnification through a filter (bandpass, 3700-5600 A) 
and onto the cathode of a photomultiplier tube. The rise-time of the photo- 
multiplier, complete with load resistor and cabling, was checked by a 
gallium phosphide photo-diode that was driven by a square-wave generator. 

The 1/e rise time was close to 1 microsecond. 

Gas Mixture 

The test mixture was prepared by the method of partial pressures in a 
stainless-steel tank. Oxygen and hydrogen were research grade gases. Carbon 
monoxide was CP grade and was cold-trapped to remove carbonyls. Argon had 
a stated purity of 99-99 percent and was cold-trapped to eliminate water 
vapor. After preparation, the mixture was allowed to stand for a week to 
insure homogeneity. 

RESULTS AND DISCUSSION 
Experimental Data 

Figure 2 shows the oscilloscope record obtained for the run at 90 torr 
initial pressure. The clean, noise-free character of this record is com- 
pletely typical of all the data. After arrival of the shock wave at the 
observation point, indicated by the deflection of the middle trace, no light 
was observed for some time. Then, a brief and almost symmetrical pulse of 
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flame -band radiation appeared, reflecting the anticipated rise and fall of 
atomic-oxygen concentration. The lower trace shows the th in-film signal 
at higher gain; the continuous rise following the initial jump shows that 
the boundary layer was thoroughly turbulent, as desired. Had there been 
any appreciable laminar portion, it would have shown up as a horizontal 
line following the initial jump**. 

As stated earlier, the two experimentally observable properties to be 
used in the evaluation of boundary- layer effects were the exponential time 
constant, eC , and the peak- time, / tT , of the atomic oxygen concentration. 
These were obtained from semi-logari thmic plots of light intensity against 
laboratory time. Figure 3 is such a plot, made from data read off the 
original of Fig. 2. The we 1 1 -defined exponential part of the record com- 
prised 1.5 to 2 decades of rising light intensity in all runs. The slope, 
jC t of a line drawn by eye through this linear portion of Fig. 3 and 
the peak-time, 'CT , also picked off visually, are given in Table I together 
with the results of the runs at the other pressures. 

Analytical Data 

In addition to experimental results, Table 1 contains the following 
calculated data. 

Shock-front temperatures were computed from the shock velocity measured 
by the two pickups straddling the window. This computation employed a 
machine program in which full thermal equilibrium was assumed, real-gas 
thermodynamic data were used, and chemical composition was assumed frozen. 
Other results of the shock calculation, while not listed, were needed as 
inputs for the analysis; these included the pressure, density, and velocity 
of the shocked gas. 
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As already described in the section on Analysis, the two methods of 
treating experimental data were applied to the 120-torr run to obtain two 
sets of rate constants that gave close fits to the observed eC and 'tT . 
These rate constants were then used to obtain the predicted constant-area 
and varying-area values of dC and T* that are listed for each run at the 
lower pressures. 

The values of l m , needed to define the apparent area change in the 
varying-area analysis, were calculated by formulas^ that apply to the case 
of all-turbulent boundary layers. This calculation required the use of a 
quantity, designated (3 0 in reference 3, which is tabulated there for argon 
and for air. Inasmuch as the mixture contained 88.5 percent argon, the 
value of (3 0 for that gas was used. 

Comparison of Analyses 

The two methods of analysis are most readily compared in terms of the 
ratios of predicted to observed values of and 'C' . Barring experimental 
errors, a completely successful analysis should produce values of these 
ratios that are as close to unity for the lower-pressure runs as they are 
for the 120-torr run. 

Figure 4 shows the results. It is immediately evident that the varying- 
area analysis based on Eq. (1) is superior to the conventional one in which 
boundary-layer effects are neglected. The constant-area method introduces 
a discrepancy with experiment that increases as pressure decreases from the 
match point at 120 torr. 
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This pressure effect can be understood by considering Fig. 5. Here 
the computed effects of the boundary layer are shown for the 90-torr run, 
which is typical of the others in this respect. At low pressure, where the 
reaction rate is reduced and the light pulse tends to be delayed and spread 
out, there is a counteracting effect due to the rising residence time and 
rate constant. Hence, and / TT do not change with pressure as mych as 
they would in the absence of flow non-uniformities. 

It should be noted that the effects plotted in Fig. 5 are gas-dynamic 
and not due to any peculiarity of the chemistry that is involved in this 
particular experiment. Virtually identical curves result if the calcula- 
tions are repeated with all rate constants set equal to zero. 

Shock-Contact Surface Separation 

The success of the vqrying-area analysis in accounting for the experi- 
mental results suggests that the essential ingredient in the analysis, 1 m , 
should be examined more closely. Although the immediate aim of this work 
was to exhibit the effects of flow non-uniformities in a kinetic experiment 
and to ascertain that Mirels 1 approach could deal with them, the ultimate 
use of the results should be to Improve the accuracy of shock-tube kinetic 
data. In this latter regard, the use of calculated V* (calculated, 
moreover, using values of for pure argon) introduces an element of 
uncertainty. It was therefore decided to measure the lengths. 

This was done by setting up an infrared monochromator with an indium 
antimonide detector coupled to an oscilloscope, to observe 4.7-micron 
emission from CO in the gas mixture. Runs were made at 50, 70, 90, and 
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120 torr, over the same small range of shock velocities as before. A sharp 
jump in emission corresponded to the arrival of the shock at the window 
and an abrupt drop to the arrival of the contact surface. The measured 
separation length, 1, was obtained by multiplying the duration of infrared 
emission by the shock velocity. 

The measured l's were within 20 percent of the Vs calculated for the 
same pressures. Thus, it is clear why the varying-area analysis succeeded 
so well. Nevertheless, the experiments were re-analyzed using l's read 
from a line faired through a plot of 1 against initial pressure (see Table 
11). The same rate constants (Eqs. 4(b) and 5(b))were used. The results, 
listed in Table II as ratios of predicted to observed values of oC and , 
show by comparison with corresponding data based on l m that use of the 
measured l's improves the agreement inoC" and leaves 'TT virtually unchanged. 

Also in Table II are results showing the consequences of applying the 
theory-' of separation-length development, according to which calculation l's 
should have been much less than l m at the window position of the shock tube. 
Re-analysis with these predicted separations required re-fitting the 120- 
torr run; this yielded a new U -2 15 percent smaller than before. The results 
at the lower pressures (Table II) are clearly inferior to those obtained 
wi th measured l's. 

Whether the unexpectedly rapid approach to l m is a general rule or 
peculiar to this tube is unclear. Until the situation is clarified, it will 
be necessary to use measured separations in the reduction of shock-tube data. 
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Implications for Shock-Tube Kinetics 

The experimental results obtained in this experiment can be viewed as 
data designed to determine the rate constant to which the results are most 
sensitive, namely, k 2 . From this point of view, t(ie effect of applying 
boundary- layer corrections is two-fold. First, the rate constant at 120 
torr is 30 percent less than the value obtained when corrections are 
neglected; and second, repeat determinations at lower pressures yield values 
with an average error of only + 4.3 percent, while the conventional treat* 
ment of the data produces a pressure-dependent rate constant. 

This suggests that much of the shock-tube kinetic data in the liter- 
ature is wrong in some degree. It is also likely that some of the scatter 
and some of the discrepancies in the results of different investigators, 
commonly noted in shock-tube work, can be attributed to the neglect of 
boundary- layer effects. Consider, for example, what would have happened if 
only the runs at 50 and 120 torr had been made and analyzed in the usual 
way. Assuming that the activation energy were known so as to take account 
of the difference in shock-front temperatures, the result would have been 
two values of the pre-exponential part of kz that differed by 25 percent. 
This would no doubt have been accepted as experimental scatter. It is 
apparent that an even larger difference could have been found if these two 
runs had been made in separate shock tubes, especially if the tubes differed 
in length and diameter. 

The size and the sign of the errors in existing data depend in a com- 
plicated way on the size of the tube, the position of the observation point. 
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the gas mixture and pressure used, and the type of measurements made, as 
well as on the rate and activation energy of the reaction. Rate constants 
derived from measurements close to the shock front are obviously the least 
suspect. Thus, the results of an investigation that covered a temperature 
range are likely to be more reliable at the high end of the range. At the 
low end, where rates are low and observations must be extended far behind 
the shock, large errors are possible, especially if the activation energy 
is large and therefore causes the rate to grow rapidly because of the 
rising temperature during the observation period. Although the error incurred 
by neglecting the flow non-uniformities in such an experiment may be either 
positive or negative, in most cases it will be positive (Ref. k describes 
a less-common case in which the error was negative). Therefore, a general 
tendency must exist in the literature for low- temperature rate constants to 
be too large, while high- temperature values are more nearly correct, so that 
activation energies will tend to be too low. An example of this behavior is 
the activation energy of (II ), which was found to be 16.3 kcal/mole with and 
only 11.9 kcal/mole without corrections for nonuniform flow^. 

CONCLUSIONS 

Changes in temperature, density, and residence time induced by turbulent 
boundary layers in shock tubes have been shown by careful experiments on a 
simple chemical system to significantly affect both the accuracy and repro- 
ducibility of kinetic results. Existing boundary- layer theory, together 
with measured separations between shock and contact surface, satisfactorily 
accounts for these effects. The results show that much of the existing data 
based on the conventional assumption of uniform flow must, in varying degrees, 


contain some error. 
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TABLE I. - OBSERVED AND PREDICTED 


\ Initial pressure, torr 


110 

Shock velocity, mm/Jusec 

1.042 

1.045 

Shock-front temperature, °K 

1128 

1133 

Limiting separation, l m , cm 

151.0 

147.7 

Exponential time constant, oCxlO""^, sec"^ 



(laboratory time scale) 



Observed 

6.77 

6.37 

Predicted, varying area 

6.66 a 

6.40 

Predicted, constant area 

6 . 58 a 

6.31 

Peak time, T, //sec (laboratory time) 



Observed 

193 

206 

Predicted, varying area 

1 99 a 

203 

Predicted, constant area 

I87 a 

194 


a Fitted by trial-and-error changes in k2 


IT 


SAL TIME CONSTANTS AND PEAK TIMES 


100 

90 

80 

1 

60 

50 

1.04J 

1.050 

1.052 

1.042 

1.045 

1.064 

wm 

1141 

1144 

1128 

1133 

1152 

144.2 

BEB 

136.4 

131.9 

127.0 


6.53 

5.67 

5.37 

4.83 

4.69 

4.65 

6.35 

6.22 

5.97 

5.17 

4.93 

5.03 

5.89 

5.63 

5.19 

4.04 

3.69 

3.54 

205 

226 

233 

270 

311 

293 

211 

219 

230 

274 

297 

289 

207 

219 

238 

306 

340 

345 















































TABLE II. - COMPARISON OF VARY! NG-AREA ANALYSES USING DIFFERENT 


SHOCK-CONTACT SURFACE SEPARATIONS 





Initial pressure, 

torr 



A v. deviation 


120 

110 

100 

90 

80 

70 

60 

50 

f rom un i ty 

Measured separation, 1 , cm. 

144.5 

145.8 

147.0 

147.5 

D 

146.6 

145.0 

142.3 


oCCcalc'd. )/cC (obs, ) 

1.003 

1.043 

.979 

1.093 

H 

1.048 

1.023 

1.022 

.043 

X (ca Ic'd.) At (obs.) 

1.016 

.986 

1.038 

.974 

H 

1.048 

.975 

1 . 02 ? 

.025 

Limiting separation, 1 m , cm. 

151.0 

147.7 

144.2 

140.5 

136.4 

131.9 

127.0 

121.3 


oC (ca Ic'd.) (obs.) 

.983 

1.003 

.973 

1.095 

1.112 

1.069 

1.052 

1.082 

.057 

X (ca Ic'd.) /X (obs . ) 

1.032 

CO 

< J \ 

1.029 

.969 

.987 

1.014 

.956 

.987 

CM 

O 

* 

Calc'd. separation at 










window, 1 w , cm. 

95.1 

94. 1 

93.0 

91.8 

90.6 

89.4 

87.5 

85.2 


©C(calc'd.)4»C (obs.) 

1.002 

1.080 

1.019 

1.158 

1.184 

1.156 

1.153 

1.170 

.115 

X (calc'd.)/^ (obs.) 

1.018 

.971 

1.010 

.947 

.961 

.981 

.900 
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Figure jl. ( - Illustration of sensitivity of light-emission profile to changes in shock 
velocity. Calculated by constant-area method for ^/(^/CO/Ar = 5/0.5/6/88.5 
mixture at 90-torr initial pressure, processed by a shxk with velxity required 
to heat the gas to 1100° K and by shxks 1 percent faster and slower. 
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Figure?;. - Experimental record obtained from H^Og/CO/Ar = 5/0.5/6/88.5 mixture 

at 90-torr initial pressure, shocked to 1141 °K. Upper trace: flame-band rad- 
iation. Center trace: thin-film record showing shock arrival at observation 
point. Lower trace: amplified thin-film record showing immediate onset of 
turbulent boundary layer. 



Figure 3. - Semi-logarithmic plot of flame-band emission shown in figure 2. 
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Figure 4(a). - Comparison of predicted and observed exponential time 
constants. O, a( predicted, varying a rea)/a(observed). •^(pre- 
dicted, constant area)Za(observed). 
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Figure 4(b). - Comparison of predicted and observed 
peak times. O, tfpredicted varying area)/x(observed). 
•, T(predicted, constant a rea/-r(observed). 
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Figure 5.. - Variation with time of post-shock properties for 90-tor r run shown in 
figure. 2. Ratios of properties calculated by varying-area analysis to those cal- 
culated by constant-area analysis. 



